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Abstract 

Type I Hermite-Pade polynomials for a set of functions /o, fi, ■ ■ ■, fs 
at infinity, Q n , o, Q n , i> • • •, Qn,s , is defined by the asymptotic condition 


Rn{z) 


(Qn,ofo+Qn,lfl+Qn,2f'2 + ' ' '+Qn,sfs)(z) ~ O 



Z —»• 00 , 


with the degree of all Q n ,k < n. We describe an approach for finding 
the asymptotic zero distribution of these polynomials as n —> oo under 
the assumption that all ffs are semiclassical, i.e. their logarithmic 
derivatives are rational functions. In this situation R n and Q„. t kfk 
satisfy the same differential equation with polynomials coefficients. 

We discuss in more detail the case when fk s are powers of the 
same function / (fk = for illustration, the simplest non trivial 
situation of s = 2 and / having two branch points is analyzed in 
depth. Under these conditions, the ratio or comparative asymptotics 
of these polynomials is also discussed. 

From methodological considerations and in order to make the sit¬ 
uation clearer, we start our exposition with the better known case of 
Pade approximants (when s = 1). 


1 Introduction and main results 

1.1 Type I Hermite-Pade polynomials 

For s € N let f = (/o, /i, • - ►, fs) be a vector of analytic functions defined 
by their Laurent expansions at infinity, 

OO r 

/fc( z ) = E fc = 0,l,...,s. 


1 



For an arbitrary n 6 N the associated vector of type I Hermite-Pade (HP) 
polynomials, Q n ,o, Qn, 1 , Qn,s, corresponding to f, is defined by the 

following asymptotic condition: 


Rn(z) (Qn,o/o+Qra,l/l+Qra,2/2 + ' ' ’PQn,sfs) — & 

where 

Qn,k 6Pn, k = 0,1,..., s, and not all Q n ^ k = 0. 



z oo, 

( 1 . 1 ) 


We denote by P n the space of algebraic polynomials with complex coefficients 
and degree < n. Function R n defined in (1.1) is called the remainder. 

In order to avoid unnecessary complications we will assume hereafter 
that all /fc(oo) ^ 0; with this condition, and since the division of (1.1) by 
fo preserves this asymptotic relation, we can also assume without loss of 
generality that fo = 1. We follow this convention in the rest of the paper. 
Nontrivial polynomials Q n) k satisfying (1.1) always exist but, in general, 
they are not defined uniquely. To avoid an essential degeneration we will 
also require that functions /i,..., f s are rationally independent. 

Construction (1.1) for the case fk(z) = e k ^ z was introduced by Hermite 
in 1858, who used it to prove that the number e is transcendental. The 
particular case s = 1 corresponds to the (diagonal) Pade approximants, and 
the general construction is called Hermite-Pade (HP) approximation to a 
collection of functions fo, ■ ■ ■ , fs- It plays an important role in Analysis 
and has significant applications in approximation theory, number theory, 
mathematical physics and other fields. For details and further references 


see mEunuMUMusaEQ]. 

A more general classes of type I HP polynomials may be introduced by 
allowing arbitrary distribution of powers of polynomials Q n k\ in this pa¬ 
per we keep all their degrees equal (the so-called diagonal case); neither we 
consider the associated construction of type II Hermite-Pade polynomials 
(see e.g. [2H Ch. 23], [38] Chapter 4], as well as [I, [5] 39] [59] for defini¬ 
tions). There is a well known formal relationship between the type I and 
type II HP polynomials [Ml [55] that can be elegantly expressed in terms 
of the associated Riemann-Hilbert problem m- However, this algebraic 
connection does not mean that an explicit relation between the asymptotic 
behavior of both classes of polynomials exists. At least, the authors ignore 
any rigorously proved result establishing this correspondence for the HP zero 
asymptotics for a sufficiently wide class of functions. It is worth pointing out 
also that the zero distribution of the type II Hermite-Pade polynomials is 
currently intensively studied by A. I. Aptekarev and coauthors (see e.g. 0)- 
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The asymptotic theory of the Pade approximants is fairly well devel¬ 
oped after the works of Stahl and Gonchar-Rakhmanov in the 80ies, see 
e.g. nzmaEiiEa. However, the analytic theory of the HP polynomials 
(in particular, the asymptotics of Q n ^ and of the remainder R n as n —> oo) 
is still in its infancy. In particular, we lack any general analogue of the 
Gonchar-Rakhmanov-Stahl (GRS) theory for HP polynomials. Some situ¬ 
ations are rather well understood, and some particular classes of functions 
f have been studied in depth. This is the case of the so-called Markov-type 


functions (that is, when the coefficients f mt k in (1.1) are moments of positive 
measures on the real axis), and of functions with a “small” set of branch 
points. In this paper we focus mainly on the second class and more partic¬ 
ularly, on functions with a finite number of branch points, see the definition 
below. 

For a polynomial Q n E P n we denote by fJ,(Q n ) its normalized zero 
counting measure 

^ n n 

KQn) = - y~] Sx k where Q n (x) = c n JJ(1 - x/x k ), c n / 0 


k =1 


k =1 


(if deg Q n < n we assume that n — deg Q n zeros of Q n are at infinity, x k = 
oo). Thus, we associate with each polynomial Q U)k (n E N, k = 0,1,..., s) 
defined by (1.1) its normalized zero counting measure 


Hn,k — /^(Qn.fc) 


( 1 . 2 ) 


and study the weak-* convergence of the sequences {/r n ,fc} as n —> oo. 

This problem can be reformulated in the following terms. For a (finite, 
Borel, and positive) measure cr on C we denote by 

c 


its Cauchy transform. Recall that the weak convergence of a sequence of 
measures fx n to a measure /r implies the convergence of to C^, both in 
the plane Lebesque measure m 2 and in L p , p > 1, on compact subsets of C. 
Since 


Cn,k ( 2 ) 


C^ k (z) 


1 Q'n,k ( z ) 

^ Qn,k ( 2 ) 


(1.3) 


we are equivalently interested in the behavior of C n k as n —> 00 . Hence, 
the proof of convergence of these sequences and the description of their 
limits constitute an extension of the theorem of Stahl to the context of HP 
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polynomials. In this generality this problem is completely open, and even 
a reasonable approach to its solution is not clear. The most general results 
so far have been obtained for the case of the type II HP polynomials for 
two functions (s = 2), each one with two branch points [3], and even these 
results are a consequence of the strong asymptotics established using the 
Riemann-Hilbert method, which is a clear overkill if we are only interested 
in the zero asymptotics. 

Thus, here we restrict our attention to a specific class of functions, gen¬ 
eral enough to be interesting, but for which we can put forward a strategy 
for studying the weak asymptotics. Namely, for a fixed set A = {ai,..., a p } 
of p > 2 distinct points let 


{ v p 1 

f(z) = f(z-,a):=Y[{ z ~ a j) aj - aj G C \ Z, ^ aj = 0 > . 

j=i j=i ) 

(L4) 

By £ we denote the union of all classes £.4 for all finite sets A. Observe 
that £ is a subclass of the so-called semiclassical family, see the definition 
below. Since each / G £ is regular at infinity, for convenience we always fix 
its branch there as /(oo) = 1 . 

With this convention the following result holds true, which shows that 
we can associate with the HP polynomials for a vector of functions from £ a 
linear differential equation with polynomial coefficients, whose fundamental 
system of solutions can be build using the Hermite-Pade polynomials: 


Theorem 1.1. For any vector f = (/o = 1, fi, ..., f s ) of rationally inde¬ 
pendent functions fk G £, k = 1 ,..., s, and any n £ N, each of the following 
s + 1 functions 

Qn, 0 ; Qn,lfli ■ ■ ■ i Qn,sfs 


defined by ( 1 . 1 ) satisfies a linear differential equation 


n s (z)w (s) + + • • • + nfizfw' + n 0 (z)w = o, (1.5) 


with polynomial coefficients nfc(z) = n k,n(z) depending on n, whose degrees 
are jointly bounded by a value depending only on the number of branch points 
of the component fk of f. 


It follows from the definition that the reminder R n defined by (1.1) is 
also a solution of the equation (1.5). 

Linear ordinary differential equations (ODE) with polynomial coeffi¬ 
cients are one of the central topics in classical analysis, and in particular, 
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the problem of existence of polynomial solutions for a given equation is well 
known. For a second order equations these polynomial solutions are called 
the Heine-Stieltjes (HS) polynomials (see [55] § 6.8]). For relatively recent 
important developments in the theory of HS polynomials see, e.g. 133. 35 TT) . 
and also [13J, 50, j56j EH 58] for the special case of the so-called Heun poly¬ 
nomials. 

From this perspective, Theorem 0 can be viewed as a construction of a 
class of higher order differential equations, which a priori have polynomial 
solutions, as well as a fundamental system of “quasi-polynomial” solutions 
of the form Q n ,kfk, for fk G £, that constitute a direct generalization of the 
HS polynomials. 

It is based on 


The proof of Theorem |1.1| is presented in Section 3.1 


identities for the Wronskians (method which goes back to Riemann and 
Darboux, see m for further references) and strongly relies on the fact that 
any function / G £ is semiclassical: it satisfies the Pearson-type equation 


a 4 


p p 

3 j =1 3 

or in other words, / is a solution of the first order differential equation 
Af — Bf = 0 where the roots of the polynomial A are the branch points of 
/ and deg-B < degH — 2. This method is well known for Pade polynomials 

mm- 

For a special choice of the vector /, which goes back again to the original 
ideas of Hermite, we can be more specific. Namely, let / G £; a vector of 
the form 

f = (i,/,/ 2 ,...,r) (i.6) 

is a particular instance of the so-called Nikishin system (of functions). Con¬ 


dition (1.1) defining the type I Hermite-Pade polynomials for f now takes 
the form 


R, 


i(z) (Qn,0 + Qn,lf + Qn,2/^ + - ‘ ' + Qn,sf b ) (z) — O T 


sn-\-s J 


OO, 


(1.7) 

with the assumptions that all deg Q n ^ < n and not all Q n> k are = 0. 

Obviously, Theorem 0 applies in this case too, so we get an ODE of 
the form (1.5) for the Nikishin system (1.6). However, even in this situation 


the number of accessory parameters in the coefficients of the equation (1.5) 


makes the problem virtually intractable. Thus, we consider in more detail 
the case s = 2, i.e., 

/ = ( 1 ,/,/ 2 )- ( 1 - 8 ) 
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Let A = {ai,..., a p }, p > 2, be the set of the pairwise distinct branch points 
of /, so that / G £_ 4. For the sake of simplicity, we will assume initially 
that deg Qn : k = n for all k G {0,1, 2}, and that the order of the zero of R n 
at infinity is exactly 2n + 2. In the established terminology this means that 
such n G N are normal indices. If / is such that every n G N is normal, we 
say that the system (1.6) is perfect (see mm)- 


Theorem 1.2. Let A = {ai, ..., a p }, p > 2, be a set of pairwise distinct 
points, f G £_4, and Q n ,k, k = 0,1, 2, be the type I Hermite-Pade polynomi¬ 
als for the system f = (1, /, f 2 ), under the assumption that these functions 
are rationally independent. 

Then functions Q n ,kf k , k = 0, 1 , 2, constitute three independent solutions 
of the differential equation 


A 2 Hw'" + A{3(A'-B)H-AH'}w" -3(n-l)(n+2)iV + 2n(n 2 -l)Gw = 0, 

(1.9) 

( 1 . 10 ) 


where 


A (z) = Jl(* “ °J')> 

3 = 1 


B = Af /f G Pp_ 2 . 


If n G N if normal, then H = H n (z ) = z 3p 6 + • • • 
2 5p-8 + ... € P 5p _ 8; and G = Gn(z ) = z 5p ~ 9 + 
polynomials dependent from n. 


£ ^3^-6; F — F n (z) — 
■ ■ ■ G P5 P -9 are some 


It follows by definition that the remainder R n is also a solution of (1.9). 


Remark 1.3. Although the electrostatic model for the zeros of Pade polyno¬ 


mials (at least, in the case of the class C) are well known (see Remark 2.3 
below), an analogue for HP polynomials, satisfying the third order ODE 


(1.9), is currently an open problem. 


Zero asymptotics of the HP polynomials is highly non-trivial (see e.g. the 
results of the numerical experiments and their discussion in j29j]). In the 
situation of Theorems 1 1.11 and II.21 its derivation is based on the substitution 
w n = exp (n f z v n (t)dt) (equivalently, v n = ^ w' n /w n ), where w n is a solution 
of the ODE. This reduces the equation to a “compact” form (generalized 
Ricatti), for which all limit equations are purely algebraic. This argument 
can be carried through along convergent subsequences of functions C n ^ (or 
equivalently, of weakly convergent subsequences of zero counting measures). 
However, for p > 2 the existence of non-trivial set of accessory parameters 


in the coefficients of the differential equation (H and G in the case of (1.9)) 
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does not allow to show that the limit exists along the whole sequence n € N, 
without appealing to some deeper arguments. 

The situation in its full generality is so complex that we start our dis¬ 
cussion in Section [ 2 ] with the case of Pade approximants (s = 1). In this 
situation the accessory parameters still exist and constitute a problem, but 
the underlying Riemann surface governing the asymptotics is hyperelliptic, 
which makes the analysis if not simple but at least tractable. 


The second case that can be fully understood is of (1.8) with p = 2, that 
is, when function / has only 2 branch points. Without loss of generality 
they can be clearly taken to be ± 1 , so that we consider 


f(z ) = f(z;a) = 


z - 1 
z + 1 


2a € C \ Z. 


( 1 . 11 ) 


Observe that / can be extended to a holomorphic function in C \ [—1,1] 
satisfying /( 00 ) = 1 . 

For this function we consider the Hermite-Pade polynomials of type I 
for the system f = ( 1 , /, / 2 ), that satisfy 


Rn(z) (Qn ,0 + Qn,lf + Qn,2f~)(z) — O ^ 


1 


2;2n+2 


OO, 


( 1 . 12 ) 


with the assumptions that all deg Q n ,k < n and not all Q n k are = 0 . 


Now the statement of Theorem 1.2 can be made more precise: 


Theorem 1.4. Let f be the holomorphic branch of the function (1.11) at 
infinity, normalized by the condition /( 00 ) = 1 , and let Q Ut k, k = 0 , 1 , 2 , be 
the type I Hermite-Pade polynomials for the system f = (1,/,/ 2 ). Then 
Q n , 0 , Qn,if and Q n ,2f 2 (and hence, also the remainder R n ) satisfy the same 
differential equation: 


(z 2 — 1 ) 2 w"' + 6 (z 2 — 1 )(z — a)w" 

- [3(n - l)(n + 2 )z 2 + 12<xs - (3n(n + 1) + 8 a 2 - 10)]to' 

+ 2 [n(n 2 — 1 )z + a(3 n(n + 1) — 8 )] w = 0. (1-13) 


Furthermore, polynomial Q n ,2 is a solution of the differential equation 

(z 2 — 1 ) 2 w'" + 6(z 2 — 1 )(z + a)w" 

— [3(n — l)(n + 2 )z 2 — 12 az — (3 n(n + 1) + 8 a 2 — 10)]w/ 

+ 2 [n(n 2 — 1 )z — a(3n(n + 1) — 8 )] vj = 0. (1-14) 
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Remark 1.5. Observe that the ODE (1.14) is obtained from equation (1.13) 
by replacing a by —a. 

The absence of accessory parameters in the coefficients of the ODE in 
Theorem 1.4 allows us to obtain an explicit formula for the limit zero dis¬ 
tribution of the HP polynomials. In order to formulate the assertion, it is 
convenient to introduce a different branch of the function / from ( 1 . 11 ); 
namely, let in what follows 


fo{z) = 


1 - 2 
1 + 2 


2 GC\F, F = R\(- 1,1), 


(1.15) 


denote the holomorphic branch in C \ F, fixed by /o(0) = 1. Function 

Pn{z) = Q n ,i(z) +2cos(air)Q n} 2(z)f 0 (z), 

defined and holomorphic in C \ F, will play a prominent role in the con¬ 
vergence of HP approximants (see Theorem 1.8 below). Meanwhile, observe 
that 

Pn(x) = Q n ,i{x) +2cos(an)Q nt 2(x) ( f + + f~) (x), x £ (-1,1), 


where / + (resp., f~) are the boundary value of the function ( 1 . 11 ) on (— 1 , 1 ) 
from the upper (resp., lower) halfplane. 

We have: 


Theorem 1.6. Under the assumptions of Theorem l.f, let a £ M, |a| £ 
(0,1/2). Then all the zeros of polynomials Q n ,k belong to M\ [—1,1], while 
p n has at least 2n + 1 zeros on (—1,1). 

Moreover, there exist two unit measures on M, A supported on E = 
[—1,1], and v supported on F = M\(—1,1), such that the normalized zero 
counting measures p n ,k, k = 0,1, 2, defined in (1.2), converge to v, while the 
normalized zero counting measure of p n , 




p„ (x)=0 

(with zeros accounted according to their multiplicity) converges to A. 
Measures u and A are absolutely continuous, with the densities 


v (x) = — 


1 


1 


1 


2 ?r y/x 2 — 1 l -^/|x| — 1 y/\x\ + 1 


X (x) = — 


+ 


47T f/l - x 2 V^f-a: s/TT 


x 


, x £ M \ (— 1 , 1 ). (1.16) 
x£(-l,l). (1.17) 
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Remark 1.7. If a. £ C, |a| £ (0,1/2), then not necessarily all zeros of Q n ,k 
are real; however, the asymptotics (1.16)—(1.17) remains valid. 


It follows from the results and methods of mm that the limit zero 
distributions A and v above can be characterized by an equilibrium problem 
involving mixed potentials, as follows. For a measure g on C we denote by 


V*(z) = J log 


its logarithmic potential. Furthermore, for F = R \ (—1,1), let gp(z,t) be 
the Green function for the domain 12 = C \ F with pole at t: for a measure 
H supported in 12 , 

Gp{z)= [ g F (z, t) d/j,(t) 

Je 

defines its Green potential (with respect to 11). 

Then A in Theorem |1.6| is the unique probability equilibrium measure 
supported on the interval E = [—1,1] for the mixed Green-logarithmic po¬ 
tential, characterized by the following identity: 


3V x (x) + Gp(x) = 7 e = const, x £ E. (1.18) 


Furthermore, v in Theorem 1.6 is the balayage of A from 12 onto F (see 
e.g. [25] or [48] for the definition of balayage). 

Equivalently, v is the probability equilibrium measure for the mixed 
Green-logarithmic potential with respect to E, but now in the external field 


i/>(x) = 3gE(x,oo) = 3 log ^|x| + \/x 2 — 1^ , x £ F, 

that is, 

3V v (x) + G v e {x) + = const, x £ F, (1.19) 

where again 

g e( z ) = J 9E(z,t) dv(t), 

and g F {z,t ) is the Green function of C \E with pole at t. The equivalence 
of such kind of problems (in a more general setting) was discussed in [5]. It 
was also explained there that these problems are not a consequence of the 
vector equilibrium problems studied in mmm for two types of systems 
of Markov functions. 

Equilibrium conditions above were used in @5[S6] to investigate the 
convergence of Hermite-Pade approximants for a set of functions 1, /i, /2 
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under the assumptions that the pair /i, /2 forms a complex Nikishin system 
for a so-called Nuttall condenser ( E,F ), where the plate E is a finite union 
of real intervals and the second plate F exhibits a symmetry known as the 
5-property. In |46| (see also |32j ) it was shown that the Cauchy transform 
h = C v is a solution of a cubic equation, so that the Nuttall condenser can 
be associated with a three sheeted Riemann surface (see also [23]). 

Finally, the problem of the strong asymptotics of the HP polynomials 
seems to be open even in the simplest cases. For Pade approximants (s = 1) 
such formulas can be fairly easily obtained from the differential equation 
(see Section [ 2 ]) using the Liouville-Green approximation [Jl]. However, 
the Liouville-Green (a.k.a. Liouville-Steklov or WKB) method is not yet 
completely developed for higher order differential equations (for the case of 
two-point Pade approximants see |22|). so it is not clear how to obtain the 


strong asymptotics of Q n ,2 using WKB even in the simplest case of (1.14) 
with all its coefficients explicit. 

In the case of (1.11)-(1.12), since for |a| G (0,1/2) the singularities 
of the function (1.11) at ±1 are integrable, polynomials Q U: 2 satify some 
(bi)orthogonality relations on M \ (—1,1) with respect to the weight ((x — 
l)/(x +1)) 0 . It allows us to establish convergence of the ratio of HP polyno¬ 


mials, which as it follows from Theorem 1.6 will occur only on the complex 
plane cut along F = R\ (—1,1). 


Theorem 1.8. Under the assumptions of Theorem 1.6, as n —> 00 , 

Qn, 1 


Qn, 2 


(z) —> -2cos(mr) f Q (z), 


Qn, 0 
Qn, 2 


/o 2 (X>, 


( 1 . 20 ) 

( 1 . 21 ) 


locally uniformly in C \F. 


Remark 1.9. Similar results for Markov functions have been proved recently 
in 


For convenience, Theorems 1.6 and 1.8 assumed that |a| G (0,1/2), 
although the crucial constraint is that neither a nor 2a are integers. Indeed, 
if a, 2a G M \ Z, but |a| 0 (0,1/2), then function /q (see the definition 
in (1.15)) is no longer integrable on the interval E. By multiplying the 


defining relation (1.12) by a suitable real polynomial of a fixed degree, say 
m, we can reduce the problem to the previously analyzed situation of an 
integrable /o, at the cost of loosing m orthogonality conditions satisfied by 
Q n ^k and R n (see Section [4]). Although the arguments that lead to the proof 
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of the limit zero distribution in Theorem 1.8 still apply in this case (so (1.16) 
and (1.17) are valid as long as a, 2a € M \ Z), a finite number of the zeros 
of the polynomials Q Uy k is now out of our control. In consequence, instead 
of a uniform convergence as in Theorem 1.8 we can claim (1.20) and (1.21) 
in capacity only. 

We can extend these conclusions even further, to the case when a, 2a 6 
C \ Z; this is subject of a manuscript in preparation. 

relies on the exis- 


The proof of Theorem 1.8, carried out in Section 4.2 


tence of a Riemann surface with the so-called Nuttall’s canonical partition 
(see [39]) - As a consequence, we conclude that for the remainder R n , de¬ 
fined by (1.12), there exists the limit g(z) = lim(l/n) log \R n (z)\, which is a 
single-valued harmonic function on this Riemann surface with the prescribed 
behavior at infinity on each sheet of 91 (what we call the real g-function of 
the Riemann surface 91): 


Definition 1.10. Let 91 be a (s +l)-sheeted compact Riemann surface, and 
7 T : 91 —> C, 7r(() = z, ( € 91, its canonical projection. For zeCwe use 
the notation 7r~ 1 (z) = The real g-function g = g<x on 

91 is defined as a function harmonic on the finite part of 91 and satisfying 
conditions 


g(C) =~ slog \z\ + 0 ( 1 ), C->°o {0) , 

5(0 = log \z\ + 0(1), j = l,...,s, 

5(C (0) ) + 5(C (1) )+ ••• + <?( C W ) = 0, 


with 

7T 1 ((X)) = (oO 1 ), . . . , 00 (S+1) }- 

Such a function exists and is unique. 

Furthermore, any function G on 91 such that ReG = g is the complex 
g-function of 91, and can be expressed as an Abelian integral of the third 
kind. 


In theory, the considerations above can be extended to the most general 
setting of HP approximants defined by (1.7) or even (1.1), with /*. £ C, by 
saying that the limit of any convergent subsequence of (1/n) log \R n (z)\ is 
the real g-function corresponding to certain (s + l)-sheeted Riemann surface 
91.5+1 associated with f. Moreover, all limits of all the Cauchy transforms 
(1.3) along the same subsequence are rational functions on the same Rie¬ 
mann surface 91 s +i. 
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Finally, if JK s+ i can be uniquely prescribed by some additional conditions 
we conclude that the whole sequence (1/n) log \R n (z)\ converges. Such a 
condition could be formulated in terms of the zero-level set { g(z ) = 0} and 
branch types of functions in f. 

A complete implementation of the plan outlined above in the most gen¬ 
eral setting is rather a program for a future. Certain assertions needed for 
formal justification of this program are not completely proved at the mo¬ 
ment, and volume restrictions prevent us from presenting even the main 
details related to general case in this paper. Instead, we discuss next a ver¬ 
sion of our approach for the case s = 1, that is, for the Pade approximants, 
situation still far from being trivial. As it was mentioned, the Riemann 
surfaces related to this case are hyperelliptic and this is an essential simpli¬ 
fication. Also, all facts related to our method for s = 1 can be obtained as 
corollaries of the Stahl’s theorem. However, we intend to extend the method 
to the study of the Hermite-Pade polynomials where we lack any analogue 
of the Stahl’s theorem. It turns out that without recurring to this theorem 
even the case of Pade polynomials presents challenges. 

Remark 1.11. These results were in part announced in [32]. We wish to 
acknowledge the useful remarks of A. B. J. Kuijlaars on the first version of 
this manuscript, as well as of the anonymous referees. 


2 Pade polynomials 

We start the discussion illustrating our approach in the simplest situation, 
when s = 1 and f = (1,/) in (1.1). This corresponds to the Pade approx¬ 
imants to the function /, or rather, to its analytic germ at infinity; recall 
that without loss of generality we assume that /(oo) = 1, so that 


f(z) = i + J2 


( 2 . 1 ) 


m =1 


It is customary to use for the polynomials defined in (1.1) the notation 
Pn = ~Qo,n and Q n = Qi, n , so that we have the conditions P n ,Qn £ Tn, 
Qn ^ 0, and 


Rn{z)= ( Qnf~Pn)(z) = 


M n 

?n +1 


1 + o[- 

z 


oo. 


( 2 . 2 ) 


Polynomials P n and Q n (not uniquely defined by (2.2)), are called the Pade 
polynomials ( Q n are Pade denominators ), and the rational function P n /Q n 
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(which is uniquely defined) is the Pade approximant to / of order n. Recall 


that an index n is called normal if degQ n = n and M n 0 in (2.2). 


Pade approximants constitute a classical method of rational approxi¬ 
mation of analytic functions. They are the best local approximations to a 
power series or, alternatively, they are convergents of a continued fraction 
(Jacobi or Chebyshev fraction) to this series. A systematic study of such 
fractions was started in the 18th century by Euler, Lagrange and others, 
although the ground of the analytic theory was laid in the 19th century in 
the works of Hermite, Chebyshev and their students and followers, such as 
Pade, Markov, Stieltjes and others (see e.g. [6j and a recent review 0). 


2.1 The differential equation 

As before, we particularize our analysis to the case of a function / £ £ 4 , 
where A = {ai,..., a p }, p > 2 , and a 3 are pairwise distinct. The following 
theorem belongs basically to Laguerre Wl although its derivation, also 
classical (see Remark 2.2 below), is different from the original one. We 


include a more or less detailed proof here mainly with illustrative purposes, 
having in mind its extension to HP polynomials in Section [3j 

Theorem 2.1 (see m)- Let function f(z) £ £. 4 , A = {ai,..., a p }, p > 2 , 
be given by (2.1), with polynomials A and B defined by (1.10). If P n , Q n are 


the associated Pade polynomials of degree n, and R n the remainder (2.2), 
then there exist polynomials H n and C n , and a constant N (in general, all 
depending on n) such that functions P n , Q n f and R n are solutions of the 
differential equation with polynomial coefficients 


AH n w" + {{A' - B)H n - AH'Jw' - NC n w = 0, 

and deg(C n ) - deg (H n ) <p- 2. 

Furthermore, if the index n is normal then 


(2.3) 


N = n(n+1), H = H n {z) 


z) = z p 2 +- • 


£ P p _ 2 , C = C n (z) = z 2p ~ 4 + 


2p—4- 


Remark 2.2. In connection with this result see mm mm and [33 § 3, 
Theorem 3.1, formula (3.1)]). According to J. Nuttall [3D], says that the 
method of proof can be traced all the way back to Riemann m- 

Since /( 00 ) = 1, substitution of aj for — otj is equivalent to dividing 
of both sides of ( 2 . 2 ) by /. It follows that polynomial Q n is a solution of 
the equation obtained from (2.3) by changing sign of all ay, which basically 
replaces B by —B. 
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Remark 2.3. In terminology of [231, Pn s are Heine-Stieltjes polynomials. 
It is well known (see [3HE3]) that ( |2.3[ ) yields an electrostatic model for 
the zeros of P n ’s: they are in equilibrium in the external field created by 
the masses fixed at the zeros of A and B, plus by a number of “ghost” or 
moving charges, corresponding to the zeros of H n . 

Proof. We fix a neighbourhood U = {z : \z\ > R} of z = oo where series 
(2.1) is convergent. If w(z) = c\(Q n f)(z)+c 2 P n (z ) is any linear combination 
of functions ( Q n f)(z ) and P n (z ) then by the general theory, the associated 
Wronskian vanishes identically, 


w 

w' 

w" 



Pn 

Qnf 

PL 

(Qnf)' 

p: 

(Qnf)” 

= 0, z G U. 

(2.4) 


W[w,P n ,Q n f](z) = 


Since W(z) is an analytic function, it must vanish identically in its whole 
domain of analyticity, C \ A. Expanding the determinant along the first 
row yields the following second order differential equation with respect to 
w, solved in particular by w = P = P n , w = Qf = Q n f and in consequence, 
by w = R = R n \ 


W(z) = W 2 {z)w"{z) + Wi{z)w'(z) + W 0 (z)w(z) = 0, 
where (we omit the explicit reference to n ) 


W 2 (z) = 

-m(z) = 

W 0 (z) = 


P P' 

Qf (Qf)' 

p P" 

Qf (Qf)” 

P i pn 

(Qf)' (Qf)” 


p P' 

R R! 

P P" 

R R" 

pi pn 

R' R” 


= °l -2 


= O 


as 


as 


oo, 


oo, 




as 


oo. 


(2.5) 

( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


The second equality in each (2.6)—(2.8) is obtained by subtracting the first 
row of the determinant from the second one, and the order of zero for Wj 
at infinity is found using (2.2). 

So far, (2.5)—(2.8) is just a set of straightforward identities. But now we 
use the semiclassical character of /: we have Af = Bf where A = z p + , 

B G P „_2 and further 


A 2 f" = (B 2 + B'A - A'B)f = Bif, where B\ G 


2p—3- 
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Taking this into account we multiply equation (2.4) by A 2 // reducing it 
to an equation with polynomial coefficients 


A 2 (z) 

m 


W(z) = U 2 {z)w"+ U 1 (z)w'+ U o (z)w = 0, zeC\A. (2.9) 


Consider for instance the coefficient II 2 ; from (2.6) we obtain 


n 2 (*) = y w 2 (z) = A 


p P' 

AQ A(Qf)'/f 


(z) = A 


P P’ 

AQ AQ’ + BQ 


(z) = A(z)H(z ), 

( 2 . 10 ) 


where H is a polynomial; by the right hand side of (2.6), H G P p - 2 - In a 
similar vein it is established that 


n i = 


A 2 \ 

y ) ^ P 2 p—4+J ; ,7 = 0,1,2. 


Furthermore, an easy consequence of (2.6) and (2.7) is that W 2 = —W\. 
From here, 


Hr (z) = - ( Y = {A '~ B ^ H ~ AH '- 


( 2 . 11 ) 


Finally, for a normal index n, (2.9) has a polynomial solution P of degree 
exactly n. Substituting P into the equation one can calculate the coefficient 
of the highest power of z (which is z n+2p ^ 4: ) and find that the leading coef¬ 
ficient of n 0 is IV = — n{n + 1). 

Theorem |2.1| is proved. □ 


Remark 2.4. Arguments presented in the proof of theorem above show that 
if we do not assume normality, then condition (2.2) should be replaced by 

R n (z) = 


Mr, 


1 + o[- 

z 


00 , 


( 2 . 12 ) 


with 0 < £ n < p — 2, so that all indices n G N are almost perfect, according 
to the terminology of mi- Indeed, if in ( |2.12[ ) we assume £ n > p — 1 then by 
the arguments above we would obtain that in (2.3), degi7 < 0, that is h = 0 
and IIi = IIo = 0, which is not possible. It follows that for n' = deg P n = 
deg Q n we have n! > n — £ n and further N = n'(n' + 1) > n(n + 1) — 2 n£ n . 


Let us hnally single out the simplest non-trivial case of Theorem 2.1 


when / G £ 4 , with A = { 01 , 02 ,< 23 }. Then the coefficients for the normal 
indices in the differential equation (|2.3[) have the form H n (z ) = z — z n , 
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C n (z) = (z — b n )(z — v n ), where z n . b n , v n are some parameters, so that 
equation (2.3) becomes 


A(z)(z-z n )w ,, +((A'-B)(z)(z-z n )-A(z))w'-n(n+l)(z-b n )(z-c n )w = 0. 

(2.13) 

This was the equation used by J. Nuttall m to derive a formula of strong 
asymptotics for Pade polynomials and their remainder in the case / £ £ 4 , 
and which we discuss next in the general setting. 


2.2 Asymptotics 

The zero asymptotics of numerators and denominators of Pade approximants 
for functions with singularities constituting a set of zero capacity was derived 
in a seminal work of H. Stahl [she], see the details below. However, the 
proof of Stahl’s theorem is based on rather sophisticated potential theoretical 
arguments, and has no simplification even for the semiclassical functions 
from £, for which we have the bonus of the Laguerre differential equation 


(2.3). It is tempting to use (2.3) to study the asymptotic behavior of the 


Pade polynomials, finding a simpler proof of Stahl’s result, at least for the 
class £. 

The first attempt in this direction for non-classical situation, although 
not successful, is due to Laguerre himself [23]. A hundred years later J. Nut- 
tall m repeated the attempt, this time successfully, even before the pub¬ 
lication of Stahl’s results. Since the Laguerre differential equation is of 
order 2, it is natural to use the Liouville-Green (Liouville-Steklov or WKB) 
method, which actually yields even the strong asymptotics of the polynomi¬ 
als (at least, away from the support of the limiting zero counting measure). 
However, in the case of functions / £ £ with p > 2, equation (2.3) has 
non-trivial accessory parameters H n and C n whose asymptotics is a priori 
not clear. Their behavior was conjectured by Nuttall first, and then rigor¬ 
ously proved for p = 3 in the pioneering work m, where the crucial step 
relied on the Stahl’s theorem just appeared in [52]! Nutt all’s method was 
extended to the case of an arbitrary number of branch points in [35]. Again, 
asymptotics of the free parameters along full sequence n is obtained using 
Stahl’s theorem. 

In the approach above the goal was the strong asymptotics, and the 
Olver’s lemma m played the key role. Now, interested in the weak asymp¬ 
totics, we proceed slightly differently: we use the Ricatti substitution in 


(2.3) and take limits along converging subsequences, resulting in an alge¬ 


braic, instead of a differential, equation. 
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Namely, substitution 


Wn ~ ex ^ ( n / v n (t)dtj or, equivalently, v n = —w' n /w n (2.14) 


in equation (2.3) and division by AH n reduces it to the Ricatti equation 

- v' n = v 2 n + s n v n + r n , (2.15) 

n 

where 

1 / A' R TJ> \ N C n , 

r n = -(2.16) 


1 fA'BH' 

Sr,. = - - 7 - TW Nn and 


n V A A H„ 


n 


2 AH' 


Equation (2.15) has, in particular, solutions v n = P' n J(nP n ), Q' n /(nQ n ) + 
f'/(nf) and R' n /nR n . 

Due to the weak compactness of the zero-counting measures involved in 
the coefficients of the equation (2.15), we have that 

v v n (z) 


lim 


= lim s n (z) = 0 m 2 -a.e. 


n n 

on compact subsets of C. Furthermore, N/v? —> 1 as n —> oo. Consequently, 

Lemma 2.5. We have 

C n 


lim vi — 


AH, 


= 0 m, 2 -a.e. 


on compact subsets of C. 

In particular, this procedure will yield an algebraic (vs. a differential), 
and more precisely, a quadratic asymptotic expression. 

The sequence of rational functions C n /AH n is compact in the m 2 con¬ 
vergence, with possible limits including identically 00 and 0 ; this happens 
if some zeros of C n or H n (or both) go to infinity. To handle such an event 
one has to use the spherical normalization of polynomials, and normalize 
the term w 2 accordingly. 

Let us select convergent subsequences 


C n —> C and H n —> H as n E A C N. 


In the generic case we have deg H = p — 2, degC = 2p — 4; otherwise, the 
degrees of C and H may be reduced. Eventually, this is not important since 
we have cancellation of possible large zeros of C n and H n . More exactly, 


17 








Lemma 2.6. We have C = VH where V(z) = z p 2 + .... 


This lemma is eventually a consequence of the fact that the only possible 


singularities of any solution of (2.3) are the zeros of A 
In 


one can find the proof of a stronger assertion: for any large 
enough n we have representation C n = H n V n where \H n — H n \ = 0(l/n) on 
compacts in C provided that polynomials are spherically normalized. 

However, as it was mentioned before, all these arguments work along 
subsequences of N, and the uniqueness of the limiting polynomial V (and 
thus, of the asymptotic zero distribution) does not follow from this analysis. 
For p > 2 (see HE®]), uniqueness has been established so far recurring to 
the Stahl’s theorem on the weak asymptotics, that we briefly outline next0 
It is worth mentioning however a new approach, discussed in | [43j , which 
uses fixed points arguments to prove uniqueness of V. 

Let us recall that for a function from the class £ Stahl’s theorem asserts 
that there exist a unique (up to subsets of capacity zero) compact set F = 
Ff C C, which is a union of analytic arcs, with the following properties: 
the complement to F is connected, / is single-valed in C \ F, the jump of / 
across any arc in F is not identically zero and, finally, the S-property holds: 


dg 

dn< 


z = 


d9 (z), z G F°, 


dn- 


where g denotes the Green function of C \ F with pole at infinity, and n± are 
two oppositely directed normals to F°, where F° is the union of open parts of 
arcs constituting F. Furthermore, the Robin (equilibrium) measure of such a 
compact F is precisely the weak-* limit of the zero counting measures g{Q n ) 
for the Pade deminators Q n . It was proved also that the sequence of Pade 
approximants it n = Pn/Qn associated with / converges to the function / in 
capacity in the complement to F: the exact rate of convergence in capacity 
was also determined. 

The proof of the existence of such a set Ff relied on its characterization 
as the set of minimal capacity: 


cap (Ff) = min{cap(I ? ) : / is holomorphic and single-valued in C \ F} . 


This extremal problem is close to the classical Chebotarev’s problem 
of minimal capacity in the class of all continua F on plane containing A. 
For finite sets A it was solved by Grotsch’s in 1930 [20j. Stahl’s result 


1 A reader interested strictly in the proof of the main results (Theorems 
skip the rest of this section and move to Section [3] 


1.1 


1.81 may 


18 









is more general even for finite sets A; it is actually a theorem from the 
geometric function theory related to a version of a general moduli problem; 
see [54]. For finite sets A a simple solution of the existence problem based 
on max-min energy problem was given in [33j; see also f33] for a study of so 
called critical measures in plane which present another generalized version 
of moduli problem. Stahl’s theorem was extended later to the case of the 
existence of an external field (or to a varying orthogonality) by Gonchar and 
Rakhmanov in cn, and to more complex equilibria in [8j [33]. It is worth 
mentioning also [7] where an analogue of Stahl’s theorem for the case of m- 
point Pade approximants was proved. In this situation the external field is 
given by unit negative charges supported at the m interpolation points, and 
as a consequence in a “generic case” the corresponding S'-curve makes an 
optimal partition of the Riemann sphere into m domains centered at these 
m interpolation points. 

The proof of the convergence assertions of Stahl’s theorem is based on the 
complex (non-hermittian) orthogonality conditions for Pade denominators 
Q n , of the form 



Qn(z)z k f(z)dz = 0, 


k = 0 , 1 ,..., n — 1 , 


where integration is taken over any system of contours separating F from 
infinity. 

We can describe the extremal set Ff in terms of trajectories of a quadratic 
differential: there exists a polynomial V(z) = Vf(z) = Y\ P l Zi( z ~ v i) such that 
the quadratic differential — (V/A) ( dz ) 2 on the extended plane C is closed: 
all its trajectories, given by 


YM 

A(z) 


(dz) 2 < 0 , 


are either closed contours or critical arcs, joining poles or zeros of AV. 

Function y/V(t)/A(t) has a holomorphic branch in II = C \ F and the 
Green function for II with pole at infinity can be written as 


g(z) = KeG(z), G(z) = f \JV(t)/A(t ) dt, (a € A), (2-17) 

J a 


where branch of the root is selected by the condition 


lim zy/V(t) / A(t) = 1. 

z —S*-oo 
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Formula (2.17) shows that function g has a harmonic continuation to the 


hyperelliptic Riemann surface 91 of the function y/V A, which is convenient 
to interpret as a two sheeted branched covering over C. In fact, it is the real 
g-function of 91, while G is the corresponding complex g - function, in the 
sense of Definition 11.101 

Now, the problem is to prove that the Riemann surface 9v is uniquely 
determined by the function /. We suggest a procedure consisting of two 
steps. The first step is the determination of the family of all Riemann 
surfaces associated with all functions / £ with a fixed set A of branch 
points. This determination is made by using special properties of the g- 
functions associated with Riemann surfaces 91 originated by compacta Ff. 

It follows from Stahl’s theorem that the family of possible Riemann 
surfaces 91 is finite. Independently of Stahl’s theorem we can assert this 
fact from the following perspective (in the situation of common position). 

Consider all hyperelliptic Riemann surfaces with 2p—2 quadratic branch 
points out of which p are fixed at zeros of A and the remaining p — 2, not 
necessarily distinct and that we denote by v\,.. . , u p _ 2 , are free; let’s 

p -2 

V(z) = II(*-*)- 

i=1 


It turns out that the requirement that the derivative G' of the complex g- 
function can have poles only at A singles out only a finite number of such 
polynomials V. Note that the formula G' = yJV/A above tells us that this 
must hold indeed for </-functions originated by Stahl’s compacta Ff. We 
assert that this property is characteristic for the S-compacta for functions 
with branch points at the roots of A. 

Next, it is clear that the extremal compact Ff is the projection of the 
zero level of g<R, or, better to say, the zero level of gy\ is the lifting of Ff 
onto 91, that is, 

F* = {(e 91: g(C) = 0} = n~ 1 (F f ). 


The selection of the unique Riemann surface associated with / £ C 4 is 
made using properties of the projection of the zero level of g<R, which is Ff. 
this set of analytic arcs on the plane has to make / single valued in the 
complementary domain and also the jump of / across any those arcs must 
be not identically zero. 

The uniqueness of such a Riemann surface implies convergence of the se¬ 
quence R' n / ( nR n ) (and in consequence, of the sequence of counting measures 
for the Pade denominators) along the whole N. 
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3 The differential equation for Hermite—Pade poly¬ 
nomials 

The considerations of Section [2] will be extended here to the case of the 
Hermite-Pade approximants. 


3.1 Proof of Theorem ll.il 


From the definition (1.1) and the properties of C we know that the remainder 
w(z) := R n (z) is a multivalued analytic function in C, with a finite number 
of branch points. Any branch of this function in the neighborhood of infinity 
U = {z E C : \z\ > R} is a linear combination with constant coefficients 
of the functions Q n ,kfki k = 0,1 ,..., s. We will drop the index n in the 
notations, and write them Qkfk■ It follows that any branch of w(z) := R(z) 
in U is a solution of the s-th order differential equation 



w 

w' 

w" 

w W 


Qo 

Qo 

Qo ■ 

• . 

W{z) = 

Qih 

(QifiY 

(Qifi)" ■ 

• • (Qi/i) w 


Qsfs 

(Qsfs)' 

(Qsfs)" ■ 

• • (Q s fsY s) 


where denotes the k- th derivative of /. Observe that this is an s-th order 

differential equation with respect to w, solved in particular by w = Qjfj and 


in consequence 

, by w 

= R = R 1 

i, so that 



w 

w' 

w" 

yj( s ) 


Qo 

Qo 

Qo ■■ 

el" 

W(z) = 

Qih 

(QifiY 

(Qif)" •• 



R 

R! 

R" 

R (s) 


Expanding these determinant along the first row yields 

S 

W(z) = J2 W j(z)w( j \z) = 0, (3.1) 

j =o 

where the coefficient W 3 , j = 0,1,..., s, is the minor obtained by deleting 
the first row and the (j + l)-th column in the determinantal expression for 
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W. For instance 


W 0 (z) = 


Qo 

Qo • 

■ • Q ( 0 S) 

(QifiY 

(Qi/i)" ■ 

• • (Qifi) (s) 

(QsfsY 

(QsfsY 1 ■ 

•• (QsfsY S) 

Qo 

Qo ■ 

• • Q ( 0 S) 

(QifiY 

(QihY' • 

• • (Qifi) (s) 

R’ 

R" 

R ( s ) 


(*) 


. s(s+3) 

(z) = O z~ 2 


oo, 


where we have used that Qj € P n and the asymptotics of fk and R at 
infinity. Operating in the same fashion, we conclude that 


/ s(s+3) . A 

Wj(z) = 0\z — +J ), z-> oo, j = 0, 


(3.2) 


The assumption that fk E C implies that there exist a pair of polynomi¬ 
als, Aj and Bj, Aj monic, and deg(Afc) — deg (BA) > 2, such that 

& (z) = *L 
fk A k ’ 

and thus, 

j Ai) _ 


A Uk = B kjfk, j = 0 , 1 , 2 ,..., 
where Bkj are algebraic polynomials. A consequence is that 

^(Qkfk) U) e P, j = o,i,..., k = o,i,...,s. 
fk 


Hence, multiplying (3.1) by 


A s 

n aw. 

fc =0 Jk 


we obtain the equivalent ODE of the form 


A 


II = J2m z W k H z ) = o, 


k =0 


k =0 


where all coefficients are polynomials. By (3.2), their degrees are jointly 
uniformly bounded. 

This proves the theorem. 
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3.2 Proof of Theorem 11.21 


Using the arguments of the proof of Theorem 1.1 we conclude that 


W(z) 


w 

w' 

w” 

w'" 

Qo 

Qo 

Qo 

Qo 

Qif 

(Qif)' 

(Qif)" 

(Qif)'" 

Q2P 

(Q2P)' 

(Q2P)" 

(Q2P)" 1 


= 0, z £ U = {z £ C : \z\ > R}, 


(3.3) 

where we again omit the subindex n. From the assumption that / £ 

A = {ai,..., Op}, p > 2, it follows that this identity can be extended to 
C \ A. Furthermore, there exist algebraic polynomials B, C and D with 
complex coefficients, such that 


Af' = Bf, A 2 f" = Cf, A 3 f”' = Df, A(z) = H(z-a j ). (3.4) 


3 = 1 


As in the proof of Theorem 0 after multiplying W by A 6 // 3 we get a 
differential equation with polynomial coefficients. Actually, it is sufficient to 
multiply W by A 5 //' 3 . Indeed, derivatives of order three are all located in 


the last column of the determinant in (3.3), and therefore, each coefficient 


in first row of decomposition contains only one such term. 
Thus, we have proved that 


A 5 

— ttW(z) = n 3 ( 2 )u;" , + n 2 ( 2 )t(; /, + ni( 2 ;)t(; / + no( 2 )rc = 0, z £ C\A, (3.5) 

P 


where IF, £ P. The generic case conditions imply that deg Qj = n, j = 0, 1, 
and R(z) = M n /z 2n+2 + ..., where M n A 0- It implies that for polynomials 
IF,- we have degllj = 5p — 9 + j, j = 0,1,2,3. Thus, those polynomials 
depend on n but their degrees are uniformly bounded. 

It is easy to see that for II 3 we have representation 


n 3<» = A 2 (z)H(z), 


where H(z ) = H n (z ) = z 3p 6 + • • • £ P 3 P - 6 - Next, using the fact coefficients 


of equation (3.5) were by first row decomposition of Wronskian we obtain 


(with arguments similar to what was done in Section 2.1) a representation 
for II 2 : 


n 2 = - 


f 

A 5 


IT 


'A 5 

P 


P rr\'^_ 
P 


A 3 


H 


= A{3(A' - B)H - AH'}- 
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Finally, taking into account the fact that equation (3.5) has a polynomial 


Q o of degree deg Q o = n, as well as another solution with the leading term 
at infinity of the form f?* (2) = 1 /z 2n+2 + ..., we arrive at two equations for 
the leading coefficients of the polynomials IIi E Ps p _8 and IIo E ^ 5 ^- 9 : 


IIi(z) = -3{n-l)(n + 2)Ul(z), 
Theorem |1.2|is proved. 


n 0 («) = 2 (n - 1 )n(n + l)IIo(z). 


3.3 Proof of Theorem 11.41 

We turn to considering the case p = 2 , a\ = — 1 , 02 = 1 and f(z) = f(z ; a ) = 
((2- l)/(z + l)) a , 2 aEC\Z. 

It is directly verified that in this particular case we have H(z) = 1 , 
112 ( 2 ) = 6 A- 2 (z)(z — a) = 6 (z 2 — 1 )(z — a). 

Next, we prove that for polynomials. II1 = ni(z;a) and IIo = IT)(2; a) 
the following equalities are valid: 

IIi(—2; -a) = IIi(2; a), —n 0 (—2; -a) = n 0 (2; a), 

IIi(2; —&) = Ili(2; a) — 361(2) + 6(2 — a)6o, 

^2(2)110(2; -a) = ^2(2)110(2; a) + 62(2) + 6(2 — 0)61(2) + IIi(2; a)6 0 , 

(3.6) 

where polynomials 60 = const, 61(2), 62(2) are defined from the following 


equations (see (3.4)): 
A 2 (f 2 y = 4 af = bof 2 , 


A 2 (f 2 r = 8a(2a-z)f 2 = b 1 (z)f 2 , 

e 2 _ U ( \ r2 ' 


A\{f-)'" = 8 o{( 2 a — 2 ) 4 a — 4A' 2 (2a — 2) — A 2 }f = b 2 (z)f~ 


We write equation (1.13) in the form 


A^w"’ + n 2 w" + niir 7 + IIou; = 0, (3-8) 

where n 2 (2) = 6(2 - a)A 2 (z), A 2 (z) = (2 2 - 1 ) and IIi E P,5p-8 = P2, n 2 G 
P5p_9 = Pi. In particular, polynomial 112(2) is known. Let IIi(2) = k 2 z 2 + 
d 2 z + C2, IIo(2) = k\z + ci; in the generic case we have kj / 0, d 2 ,ci,c 2 , 
functions Q n ,o,Qn,if, Qn, 2 f 2 and R n are solutions of this equation, and 
their monic versions satisfy Q* n 0(2) = z n + ..., Q* n 2 {z) = z n + ..., R* n (z) = 
l/z 2n+2 + .... Substitute w = Qo in ( 3 . 8 ); then the left hand side of 


the equation is a polynomial. Equating its coefficient at z n+l (the leading 
coefficient) to zero we obtain 


n(n — l)(n + 4) + nk 2 + k± = 0. 


(3.9) 
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Similarly, substituting w = R* n and equating to zero the coefficient at 
1 /z 2n+1 we get another equation for £q, & 2 : 


- 4(n + l)(2n + 3)(n - 1) - 2(n + 1 )k 2 + h = 0. 


(3.10) 


Solving system (3.9) and (3.10) for ki and we obtain 

k\ = 2 (n — l)n(n + 1) = 2(n 2 — l)n, k 2 = —3(n — l)(n + 2). 


From normalization /(oo;a) = 1 and identity f(z;—a) = 1 /f(z;a), 
by dividing both sides of the definition (1.7) by / 2 we come to Q^ 2 = 

Q n 0 “’ a ^. Therefore, the polynomial Q n , 2 satisfies a differential equation 
which is obtained from (3.8), by replacing a by —a. The modified equation 
is 

A^w'" + U 2 w" + IIi w' + IIou; = 0, (3-11) 

where 112 ( 2 :) = 6 (z + a)A 2 (z), klj(z) = II j(z; —a ), j = 1 , 2 . 

Next, substitute the solution w = Q n , 2 f 2 where f 2 = f 2 E 
in (3.8) and regroup terms so that equation takes form (3.11) with solution 
w = Q n , 2 - This yields the following relation for coefficients II j and Ilj of 
equations (3.8) and (3.11): 


3 A 2 B 0 + II2 — IIt, 

3 B± + 2II2-T- + IIi = III) 

A 2 

B2 B\ Bn ~ 

~ a —^ —I - no = no, 

A 2 -A 2 A-2 

where 112(2;) = 6 (z — a)A 2 (z), 

A 2 f' 2 = 4a/ 2 = S 0 / 2 , A\f'i = 8a(2a - z)f 2 = 5i/ 2 , 
Alff = f 2 8a{(2a - z)4a - 4B 2 A' 2 (2a - z) - 7l 2 } = B 2 f 2 . 

It follows from (3.12) and ( 3.14| ) 

24 az + IIi(z) = IIi(z) = IIi(z; —a). 


(3.12) 

(3.13) 


(3.14) 


(3.15) 


Using the symmetry f{—z ; —a) = f(z; a) and the differential equations (3.8) 
and (3.11), we get 

Hi(— z\ — a) = IIi(z; a), IIo(— z; — a) = — IIo(z; a). (3.16) 
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Therefore, from (3.15) and (3.16) we obtain 


IIi( 2 ; ex) = — 12 az + C 2 , where k^ = —3(n — l)(n + 2), 

Hq{z\ a) = k\z + c\, where k\ = 2{n — l)n(n + 1). 


Taking the limit as 2 —> 00 in (3.13) and with the equation ci(—a) = — ci(a) 


(following from (3.16)) we find that 

ci (a) = 2a(3ro(n + 1) — 8). 


(3.17) 


Similarly, evaluating (3.13) at z = 0 and combining it with (3.17) it follows 
that 

€ 2 ( 0 ) = 3 n(n + 1) + 8a 2 — 10. 


Finally, from (3.6D(3.7) one can find polynomials ITi and IIo explicitly: 


IIi (z\ ol) = —3 (n — l)(n + 2 )z 2 — 12 az + [3 n(n + 1 ) + 8a 2 — 10], 
IIo( 2 :; a ) = 2n(n 2 — 1 )z + 2a [3n(n + 1 ) — 8]. 

This concludes the proof. 

4 Asymptotics of Hermite—Pade polynomials 

4.1 Proof of Theorem 11.61 


Denote E = [—1,1]. Definition (1.7) (with s = 2) yields immediately the 
following orthogonality relation 


jf (Qn,lf + Qn,2f 2 ){Qq(() d( = 0, 


(4.1) 


where q 6 P 2 ?i is arbitrary, and we integrate along an arbitrary closed con¬ 
tour encircling and sufficiently close to E. Using the construction of Pade 


polynomials (see e.g. (2.2)) we find P n .o, Pn,\ € P n , P ni i ^ 0, such that 

1 


{Pn ,0 + Pn,lf)(z)=O 


?n +1 


OO. 


From (4.2) it follows immediately that 


j> Pn,l{Qf{Qp{Q d( = 0, V P GP n _! 


(4.2) 


(4.3) 
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Since f(x) > 0 when x > 1, we have 

/±( x ) = e ±i ™/„(*), *€(-1,1), (4.4) 

where / + (resp., f ~) are the boundary value of / on (—1,1) from the upper 
(resp., lower) halfplane, and /o was defined in (1.15). In consequence, 

1 — x 


A/(x) := (/ -/ )(*) = 2*sin(a)/ 0 (*) = 2zsin(a) 


1 + x 


, *€(- 1 , 1 ). 
(4.5) 


Taking into account that |a| € (0,1/2), (4.4) can be rewritten as 

J P n ,i(x)p(x)Af(x) dx = 0, Vp G P n — 1 , 

and we conclude that polynomials P n .\ = P n coincide (up to normalization) 
with the Jacobi polynomials Pn°’ orthogonal on the segment E with 
respect to the positive weight function /q: 

f P n (x)p(x) fo(x) dx = 0, p G P„- 1 . 

JE 

In order to simplify notation, let us define also 

f{z) :=2cos(a7r)/ 0 (z), z € £1 := C \ F, F:=M\(-1,1), 


so that 


/(x) = (/+ + / )(x), * G (—1,1). 


The integrability of /q at the end points of E plus the boundary condi¬ 
tions (4.4) allow us to rewrite the relation (4.1) as 


(Qn, 1 + Qn ,2 f)(x)q(x)fo(x) dx = 0, Vq G P 2 n- 


(4.6) 


Since Q n ,k have real coefficients and both functions / and /o are positive in 
the interval (—1,1), it follows from ( |4.6[ ) that the form 

Pn(z) ■= (Qn, 1 + Qn,2f)(z) (4.7) 


has at least 2n + 1 zeros in the interval (—1,1), that we denote by x n j. Let 

2n+l 

2n+l(z ) •— j | (z X n j ). 

3 = 1 
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Since for each p £ P n -i function p n p/u 2n+ i is holomorphic in the domain 
12, by the Cauchy formula we obtain 


Qn,l + Qnfif / % = J_ I (Qn,l + Qnpfm d( = J_ f (Q n , 2 /)(C) <*C 
CJ 2 n+l 2TTZJ^ W 2n +l(C)(C - ^) 2vTZ CU 2 n+l (C) (C ~ z) ’ 

(4.8) 


where 7 is an arbitrary simple analytic contour that is contained in the 
domain 0 and that contains all the points x n j, j = 1 ,..., 2 n + 1 , and z 
inside. Also 


Q 1 / (Qn,i + Qn, 2 /)(CMC)dC 1 / (Q^/XCMC)^ (4 _ g) 

2 ni J 7 u; 2 n +i(C) 2 irij 7 u 2 n +i(C) 

It is easy to see that we can transform the contour 7 in such a way that the 
relation (4.9) takes the form 

0 = f A7 to ) = /+(„)- /-to), P 6 p„_„ 

Jf W2n+l(y) 

(4.10) 

Since 


f^iy) = e =FQ,r *2cos(a7r) 


1 - 2 / 


1 + 2 / 


, y £ F° := [— 00 , —1) U (1, 00 ], 


we obtain that A f(y) = —4?sin( 2 a 7 r)|(y — 1 )/{y + 1)|", y £ F°. Thus the 


relation (4.10) may be written in an equivalent form as 


0 = / Qn,i{y)p{y)h r 
Jf 


*(y)dy, pe 


n— 15 


(4.11) 


where the weight function h n>a (y) = |(y-l)/(y+l)|“ /u 2n+ i{y). Since u 2n +i 
is of degree 2 n + 1 , function h n ^ a is negative on (— 00 , — 1 ) and positive on 
( 1 , 00 ), and it follows immediately from (4.11 ) that real polynomial Q n p has 
at least n — 1 zeros on F°. 

If deg Q U)2 = n — 1, it shows that all its zeros are contained in the F°. 

Let deg Q n p = n. Since Q n , 2 is a real polynomial with at least n — 1 real 
zeros, we conclude that all zeros of Q n ,2 are real, and at most one of them in 
E. Assuming the existence of one zero in E leads into contradiction with the 
orthogonality relation (4.11). Thus, once again we conclude that all zeros 
of Q n p are in F°. 

Since same arguments can be applied to the function (Q n ,i/f+Qn, 2 )p/u 2n +i, 
it immediately follows that we arrive at the same conclusion about the zeros 
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of the polynomial Q U) i. Finally, we can divide both sides of (1.12) by f 2 to 
obtain that Q n fl{z] f ) = Q n ,2 A; f~ 2 )- 

Hence, we have proved that all zeros of polynomials Q n ,(h Qn.i, and Q n ,2 
are contained in F°. 

Let us remark also that by (4.8), (4.9) and (4.10), 

^ 7m ^ ^ 1 I QnAC)f(C) ^ 

(Qn,l + Qn, 2 f)(z) - U 2 n+l(z) 27T . j _ z) <K 


U 2 n+l{z) 1 f Ql, 2 (C) fit) 


Qn,2(z ) 2 iri J U2n+l((){( ~ z) 


d( 


UJ2n+l(z) f Ql,2(y) A f(y) 


Qn, 2 (z) J F U 2 n+l{y){y ~ z) 


dy , z E Cl. (4.12) 


Now we turn to the asymptotics. Divide the equation (1.13) through by 
n 3 w and rewrite it in terms of h n = w'/(nw)\ 


(* 2 - D : 


( - 3 K 1 + —h n h! n + h 3 n \ + - {z 2 — 1 )(z - a) ( h' + h 2 \ 

\ rr n J n \n ) 


3(n — l)(n + 2) 2 12 az 3n(n + 1) + 8a 2 — 10 

O ^ > Q O 


71 


n* 


n* 


+ 2 


n(n 2 — 1 )z + a(3n(n + 1) — 8) 


= 0. 


r 


Since w = Q n ,0 is a solution of (1.13), the corresponding h n takes the form 

dfj-nit) 


h n (z) = C^(z) = J 


z — t 


where y n = y(Q n , 0 ) is the normalized zero-counting measure of Q n fl. 

Now we take limit along any convergent subsequenccj^J with y n A v, so 
that h n —> h, z ^ M, and get that h satisfies the cubic equation 


(z 2 — 1 ) 2 h 3 (z) — 3 (z 2 — 1 )h(z) + 2z = 0. 


(4.13) 


We apply the standard Cardano formula to find its solution. This yields the 
following result: 

Lemma 4.1. The general solution of the equation 

(z 2 - 1 ) 2 y 3 (z) - 3(z 2 - 1 )y(z) + 2z = 0 (4.14) 

2 At this point we should assume it convergent in the extended complex plane CU {00}. 
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for z 6 C \ {— 1,1} is 


2 / 0 ) = 


1 


Y(z) + 


z - 1 Y(z) ’ 


where 


Y(z) = 


1 + z 
1 - 2 


1/3 


(4.15) 


(4.16) 


In a neighborhood of infinity there are three holomorphic and linearly inde¬ 
pendent solutions yj, j = 1,2,3, that can be enumerated in such a way that 
they satisfy 


2/iO) + 2/2 0) + 2/3 0) = 0, 


1 


1 


9l(3) = ! + 0 U 

92(2) = z +° (4 


00 . 


(4.17) 


y 3 (z) = -- + 0[ 1 


Proof. A direct substitution of the right hand side in (4.15) into (4.14) 
yields that y(z) given by (4.15) is indeed a solution of this cubic equation, 
regardless the branch of Y considered. 

Furthermore, let yi(z) and y 2 (z), z ^ ±1, be two values of (4.15) corre¬ 
sponding to two different selections of the branch of Y. Then, without loss 
of generality, we can take the value of Y(z) in such a way that 


2/i0) = 


1 


Y(z) + 


1 


1 


p 2m/3 p —2iri/3 i 

V2{z) = —Y{z) + 


z + 1 ' v_ ' 1 z- 1 Y(z) } z + 1 " v ~' 1 Z -1 Y(z)' 

Assuming that y±(z) = y 2 (z), straightforward calculations yield us into a 
contradiction. Finally, relations (4.17) are obtained by totally standard 
arguments. The lemma is proved. 


□ 


Observe that the remainder function R n , defined in (1.12), is is holomor¬ 
phic in D := C \ E and is a solution of the differential equation (1.13); in 
consequence, its analytic jump w n = A R n on E° = (—1,1) is also a solution 
of the same differential equation. But 

w n (x) = p n (x)Af(x), x E (—1,1), 


where p n was defined in (4.7) and A/ is given by (4.5). Standard arguments 
show that for z ^ F, 


Ql,2(y) A f{y) d y\ \ f Ql,2(y) A f(y) d y\ 


-1 


F U> 2 n+l(y){y - z) 2 j \ J F u 2n+ i(y)(y - z) } 


0, 


n 


00 , 
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and we conclude from ( 4. 12| ) that 

1 K _ I Pn , 1 (A/)' _ 1 4„+l , , 1 Q'n.2 


I b _ » lb j 

n/)„ n A/ 


z - 


(z) + o( 1). 


^ ^2n+l Qn ,2 

Using the result of Theorem |1.6[ we conclude that 

lim-^ = 2<7 A (z) - <T(z), zG C\R. 
n nw n 

From these considerations it follows that C v and 2C X — C u are two, 
clearly linearly independent solutions of the algebraic equation (4.14). We 
already know that 


1 


C"(z) = yi (z) = - + 0 ^ , 


1 


oo. 


Taking into account (4.17) we conclude that 


1 


1 


(2 C*-C"){z) = y 2 (z) = - + 0 ) , C\z) = -~y 3 (z) = -+O^J, 

which proves in particular that A is a probability measure supported on E. 

Since we know the support of both A and is, we can recover their densities 
using the Sokhotskii-Plemelj formulas, 

v'{z) = (z/i+O) - yi-{x)), z e supp(zz) = F = R\ (-1,1), (4.18) 

A'(z) = (y 3 +(x) - y 3 -(x)) , z £ supp(i/) = E = [-1,1]. (4.19) 


Let Y now denotes the holomorphic branch of (4.16) in C \F, determined 
by T(0) = 1. With this convention, 

lim Y(ix) = e ±7n//3 . 

X —^zboo 

Since y±(x) = C v (x) and supp(z/) = M\ (—1,1), it follows that yi(x) £ M 
for x £ ( — 1,1). This allows us to single out the expression for y\. Expanding 
at infinity we obtain the expression for all three solutions in C \ M: 


1 


y 1 (z) = C'(z) = -—Y(z) + 


1 


1 


z + 1 


z - 1 Y(z ) ’ 


y 2 (z) = (2C x - C u )(z) = -yi(z) - y 3 (z ) 
y 3 (z) = -2 C\z) = 


^2m/3 e ~2iri/3 ^ 

Y(z) + 


z + 1 


g—27ri/3 

Z + 1 


-Y{z) + 


z - 1 Y (z) ’ 

g 27ri/3 ^ 


z - 1 Y(z)' 


Irn z > 0, 


Irn z < 0. 


(4.20) 

(4.21) 

(4.22) 
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Let Y + (reps., YJ) denote the boundary values of the selected branch of Y 
on (l,+oo) from the upper (resp., lower) half plan^] Then 

Y± ( x )=e ± ^{f^ x>l, 

V x — 1 

where we take the positive values of yf, so that 


A Y(x) = Y+(x) - Y(x) = 

V x — 1 


A y^' 2) y. (x) y_(x) 


x — 1 

X + 1 


Using (|4.18p and (|4.20|) we see that 

1 


1 


Az) = -5Z7 


27TZ V X + 1 


- M 


AY (®) + 

X + 1 1 , lx — 1 


27t I r + 1 V s - 1 x — 1 V x T 1 J 


which establishes (1.16). 


On the other hand, for x G (—1,1), by ( |4.19| ) and ( |4.22[ ), 

1 


Az) = 


sin(27r/3) / 1 


2vr 


x + 1 


y(x) - 


i 


x — i y(x) 


i 


1 + X 1 

+ 


1 — X 


V3 

47r\x + lVl — x ' 1 — x V 1 + £ 


where we again take the positive values of yf-, and (1.17) follows. 


4.2 Proof of Theorem 11.81 


By (4.12), 


Ql,2(y) A f{y) d y 

Qn, 2 V ~' Ql, 2 ( z ) Jf U 2n +l{y)(y - z) ’ 


Qn, 1 / \ , r/ \ Pn ( \ w 2n+l(~ 

rr-(z) + f(z) = — \z) = 

On,2 


with p n defined in (4.7). From Theorem 1.6 it follows that 

1 Qn,2 


lPn 

np n 


Qn,2 


2 (C*-C v )(z), z€ <c\: 


3 Calculations for (—oo, — 1 ) are similar. 


(4.23) 


(4.24) 
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In the previous section we have established that 

yi {z) = Cr(z), y 2 (z) = (2C x -C I ')(z), y 3 (z) = -2C\z) 


are three independent holomorphic solutions of the algebraic equation (4.14) 
in C \ R. With this notation we conclude from (4.24) that 

1 


log 

^ tyn,2 


j teCC)- yi (OHO zO 


or equivalently, 


Pn 


Qn, 2 


(0 


1/n 


exp 


! [ ReJ Z (y 2 -y 1 )(C)d^ (4.25) 


From (4.23) and (4.25) it follows that ( 1.20|) is established if we prove that 


Re / ( 2/2 - yi)(t)d( < 0, z<£\ 


(4.26) 


Let us consider the three sheeted Riemann surface 9 I 3 of genus 0 given 
by the equation w 3 = (z— \)/{z+l). It can be realized as shown on Figure[l] 



Figure 1: Sheet structure of the Riemann surface 9 I 3 . 


A general approach for establishing (4.26) is through the analysis of 
the global structure of the critical trajectories of the quadratic differential 
5 '( z)dz 2 on IR 3 , where 


S(0 


— ( 2/2 — J/ 3) 2 (z), 

if ze ^4 1} 

- ( 2/1 - 2/s) 2 (z), 

if z € 9^ 2) 

1 

<cT 

l 

to 

to 

PnP 

if z € 
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is meromorphic on 9 ^ 3 . Here again we use the notation n : 9 I 3 —> C for the 
canonical projection of 9^3, and = {z^\z^ 2 \z^} for z 6 C, with 

z® G 9f^. 

However, in the particular case under consideration we can use a result of 
J. Nuttall; following Nutt all’s approach, it is sufficient to show that the three 
functions yj(z) give the canonical (in Nuttall’s sense, see [32 sec. 4.3.4], [461 
§6]) partition of 9^3 into three sheets in the following way: 


Re 


/ z rz rz 

y 3 (z) dz < Re / y 2 (z) dz < Re / yi(z) dz, 


(4.27) 


Alternatively, we need an abelian integral </> on 94s such that 

• Re(</>) is single valued on 9^3 and harmonic on 9 I 3 \ u| =1 oo^; 

• exp {(j>(z)) is meromorphic on 9 I 3 , whose divisor on 9 I 3 is 2oo®—oo^ — 
00 


With such a q i>, Nuttall’s canonical partition (4.27) of 9 I 3 into three sheets 
is given by 

Re<^(z^ 3 ' ) ) < ~Re<f>(z^) < Re 


The key fact we can exploit is that the algebraic curve defined by (4.14) 


is independent of a, so we can set a = 1/3 and use the Nuttall’s example 
Section 4.3.4], where he considered the Riemann surface 9 I 3 of the equation 
zw 3 = z — 1, along with the Hermite-Pade approximants to the system 
f = (1,/,/ 2 ) with f(z) = w. Nuttall showed that exp(</>(z)) = z( 1 — w) 3 , 
as well as proved thal|^] 


Recf)(z^) — Re</>(z < - 1 ' ) ) < 0, z£ 


(4.28) 


Clearly, the Riemann surfaces 9 I 3 and 9 13 are isomorphic in such a way 
that (4.28) implies (4.26). 

Finally, recall that identity f(z\—a) = 1 /f{z\a) implies, by dividing 
both sides of the definition ( |1.7[ ) by / 2 , that Q^.\ Hence, 

( 1 . 20 ) applied to f(z ; —a) yields 

77 ^( 2 ) —» —2 cos(a 7 r)/ fo(z), 


which allows to conclude (1.21). The theorem is proved. 


4 The numeration of the sheets of 9I3 in [ 39 ] is slightly different. 
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